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betcer  the  correlation.  3)  The  relationship  between  time  lag  and 
wavelength  for  possibly-connected  outbursts  appears  to  be  characteris¬ 
tically  different  for  different  objects.  4)  No  pronounced  3— am  outbursts 
on  a  tine  scale  of  one  to  a  few  days  were  observed,  but  three  sharp  drops 
or  "quenchings*  were.  A  list  of  other  quenchings  foupd  in  the  literature 
Is  included.  5)  The  amplitudes  of  outbursts  vary  as  ,  where  <k>  *" 
from  3.3  to  111  m,  a  factor  of  34.  Results  4  and  5  may  place 
significant  new  constraints  on  models. 
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I.  INTRODUCTION 


The  first  observations  of  extragalactic  sources  at  3  am  were  made  with 
the  Aerospace  4.6-m  telescope  in  1965.  At  that  tine  longer  wavelength  data 
were  revealing  that  many  sources  were  variable  and  that  the  variation  time 
scales  were  often  shorter  at  shorter  wavelengths.  The  Aerospace  Corporation 
had  the  only  facility  available  on  a  regular  and  long-term  basis  for 
monitoring  at  short  millimeter  wavelengths.  Hence  we  commenced  a  program  of 
extensive-  monitoring.  The  objective  was  to  provide  information  which  could 
add  to  the  physical  understanding  of  the  sources.  Results  through  1970  were 
given  by  Fogarty  et  al. *  This  report,  which  is  a  sequel  to  and  supersedes 
that  paper,  presents  the  results  for  3C  84,  3C  120,  3C  273,  OJ  287,  and  BL 
Lac,  the  sources  we  observed  most  frequently  from  1965  through  mid-1974. 

Since  then  the  Aerospace  telescope  has  been  utilised  for  continuum  studies  of 
planets  and  millimeter-wave  spectroscopy.  We  plan  to  resume  source  monitoring 
upon  Installation  of  a  cooled  receiver  now  under  construction. 


II.  OBSERVATIONS 


Fogarty  et  al.  *  and  Epstein  at  al.2  have  described  the  telescope  and 
receiver  and  the  calibration  and  observing  procedures.  Here  we  give  only  a 
sunary  plus  some  additional  relevant  details. 

A.  INSTRUMENTATION  AND  TECHNIQUE 

The  observations  were  aade  at  3.3  on  (90  GHz)  using  the  "on-on**  dual-beam' 
technique.  Since  March  1969  the  telescope  tracked  a  source  such  that  the 
antenna  traced  out  the  identical  hour  angle  arc  during  the  two  halves  of  the 
"on-on"  pair,  resulting  in  cancellation  of  antenna  sidelobe  and  back lobe 
pi ckup  of  the  atmosphere  and  ground.  Data  prior  to  March  1969  were  corrected 
for  incomplete  cancellation.*  Each  observation  consisted  of  «  10  hours  of 
integration  separated  into  *  2 5-mi n  segments  which  were  reduced  individually 
and  then  averaged  together;  the  Integration  time  was  generally  spread  out  over 
a  few  days  in  an  attempt  to  randomize  the  effects  of  any  unknown  systematic 
errors  which  might  vary  from  day  to  day. 

Correction  factors  to  the  measured  flux  densities  for  pointing  errors 
(not  corrected  for  while  observing)  and  antenna  jitter  ranged  from  *  1.01 
to  »  1.10.  Antenna  pointing  information  prior  to  1971  was  derived  from 
frequent  measurements  of  the  limbs  of  the  Sun;  subsequently  both  solar  limb 
measurements  and  planet  measurements  by  a  five-point-grid  procedure  (FIVE) 
were  used.  The  FIVE  procedure  is  the  preferred  method  for  observing  stronger 
sources,  such  as  the  planets,  because  the  Information  for  correcting  the 
measured  value  for  painting  errors  is  contained  in  the  measurement  itself. 

B.  CALIBRATION 

The  variation  of  antenna  gain  as  a  function  of  zenith  angle  was 
determined  from  862  FIVE  measurements  of  Venus  and  Jupiter  taken  from  1971 
through  1977.  The  derived  gain  relative  to  the  zenith  value  (e.g. ,  0.923  ± 
0.010  at  z  »  60  deg)  agrees  to  within  1  -  2Z  with  that  predicted  from 
reflector  surface  accuracy  measurements  at  various  zenith  distances. 
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The  Sun  ws  the  primary  standard  of  antenna  teaperature.  A  gas  tube  was 
the  Intermediate  standard  prior  to  Noveaber  1968,  and  a  hot  resistive  load 
thereafter;  the  aeasureaents  of  the  load  were  corrected  for  aabient 
teaperature  effects  subsequent  to  data  acquisition.  The  intermediate  standard 
was  calibrated  against  the  Sun's  90-GHr  brightness  teaperature  of  7820  by 
the  average  Intercept  method.^ 

The  relationship  of  flux  density  to  antenna  teaperature  was  determined  by 
comparing  the  average  of  29  observations  of  DR  21  taken  from  1967  to  1975 
(Fig.  1)  with  the  calibrated  90-GHz  DR  21  flux  density  of  17.7  i  0.5  Jy; 

Ullch*  has  subsequently  revised  this  value  to  17.3  ±  0.6  Jy  but  we  have  not 
asde  a  corresponding  revision  in  our  calibration.  A  factor  of  1.016  was  used 
to  correct  the  flux  density  of  DR  21  for  the  effects  of  bean  broadening  by  the 
Aerospace  telescope.  The  net  overall  absolute  calibration  uncertainty  in  our 
results  is  estiasted  to  be  •  71. 

C.  SYSTEM  CHECKS 

Because  the  observations  reported  here  are  measurements  relative  to  an 
Internal  intermediate  standard,  rather  than  aeasureaents  relative  to  standard 
sources  on  the  sky,  we  performed  the  following  checks  on  system  calibration 
and  performance: 

(I)  The  histogram  of  the  ■  25  antenna  teaperature  measurements  com¬ 
prising  each  observation  was  Inspected  to  verify  that  the  dispersion  in  the 
measurements  was  consistent  with  receiver  noise. 

(II)  Over  the  years  we  have  made  hundreds  of  observations  of  ’’blank’*  sky 
(e.g. ,  see  Fig.  8  of  Ref.  4  and  Fig.  1  of  Ref.  5)  to  confirm  repeatedly  that 
our  systea  Introduces  no  measurable  biases. 

(iii)  We  verified  the  reliability  of  Sun-derived  pointing  information  by 
comparing  planetary  brightness  temperatures  obtained  by  using  pointing 
corrections  determined  only  from  Sun  measurements  with  temperatures  obtained 
by  using  the  preferred  FIVE  procedure.  Any  systematic  differences  between 
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temperatures  obtained  before  and  after  the  introduction  of  the  FIVE  procedure 
are  ~  3X  (see  Fig.  1  of  Ref.  2). 

(iv)  The  repeatability  of  the  measurements  is  shown  in  the  display  of 
the  DR  21  data  (Fig.  1).  Even  though  these  data  suggest  an  outburst  of  2  -  3 
Jy  amplitude  during  the  second  half  of  1973,  the  scatter  about  the  mean  value 
is  entirely  consistent  with  the  l-o  error  bars;  furthermore  there  were  no 
comparable  percentage  variations  in  our  1973  Jupiter  and  Saturn  data  (see 
Fig.  1  of  Ref.  2),  Indicating  that  there  was  not  a  multiplicative  error  in  all 
of  our  data. 

(v)  He  inspected  all  the  data  for  coincidental  rises  and  declines  with 
the  aid  of  a  procedure  similar  to  that  of  Kesteven  et  al . ^  During  only  one 
interval  (approximately  S  weeks  in  the  spring  of  1972)  did  all  five  sources 
exhibit  the  same  pattern — simultaneous  rises  by  an  average  of  12Z.  These 
simultaneous  rises  are  so  unlikely  that  we  assumed  the  presence  of  an  unknown 
systematic  error  and  reduced  all  readings  during  that  Interval  by  12Z. 

(vi)  Numerous  comparisons  with  Independent  3-mm  measurements  obtained 
with  the  NRAO  11-m  telescope  by  Hobbs  and  Dent7  are  presented  in  the  lower 
portions  of  Fig.  2,  where  the  NRAO  3-nm  measurements  are  shown  as  red 
triangles  (the  height  of  each  triangle  corresponds  to  twice  the  l-o  error). 

The  Hobbs  and  Dent  data  have  been  multiplied  by  the  factor  1.047  to  account 
for  their  adoption  of  a  DR  21  flux  density  of  16.9  Jy  instead  of  the  17.7  Jy 
used  here.  The  overall  agreement  between  the  two  completely  independent  sets 
of  data  is  excellent. 
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Fig.  2.  flower  halves)  Flux  density  at  3.3  mm  vs  time.  The  error  bars  indicate  the  la  relative  uncertainties  only;  the 
estimated  overall  absolute  calibration  uncertainty  is  1%.  The  red  triangles  represent  3-mm  data  obtained  by  Hobbs 
and  Dent  (1977)  with  the  NRAO 1 1-m  antenna;  their  height  corresponds  to  twice  the  la  relative  measuring  error, 
(upper  halves)  Centimeter-wave  flux  densities  vs  time  plus  hand-drawn  representations  of  the  3. 3-mm  data.  The 
error  bars  indicate  la  relative  uncertainties.  To  avoid  overcrowding  error  bars  are  omitted  when  they  are  leas  than 
twice  the  height  of  the  symbol  (thus  the  absence  of  an  error  bar  does  not  automatically  mean  the  error  is  the  size  of 
the  symbol  or  smaller).  Note  that  the  intensity  scales  for  the  upper  and  lower  halves  are  different. 

T 


FLUX  DENSITY  (Jy)  FLUX  DENSITY  (Jy) 


60 

54 

48 

42 

36 

30 

24 

18 

12 

6 

0 


i — i — i — i — | — i — i — i — i — j — i — i — i — i — | — i — i — i — i — | — i — i — r 


0J  287 
3.3  mm 


1969  1970  1971  1972  1973 


Fio.  2.  (continued) 


1972 


1973 


1974 


1975 


FLUX  DENSITY  (Jy)  FLUX  DENSITY  (Jy) 


FLUX  DENSITY  (Jy)  FLUX  DENSITY  (Jy) 


i 


III.  RESULTS 


The  adopted  calibration  factors,  attenuation  correction  schemes,  and  data 
reduction  procedures  changed  several  tines  over  the  years.  To  assure  uniform¬ 
ity  of  data  handling,  we  re-reduced  all  9  years  of  data  for  3C  84,  3C  120,  3C 
273,  OJ  287,  and  BL  Lac,  the  five  aost  frequently  observed  sources;  the 
results  here  supersede  those  published  earlier.*  (The  data  for  many  other 
less  frequently  observed  sources*  have  not  been  re-reduced.)  The  results  are 
listed  in  Table  I  and  are  shown  in  the  lower  halves  of  Fig.  2.  The  l-o  error 
bars  shown  indicate  only  the  relative  errors  of  measurement.  The  pronounced 
reduction  in  the  errors  beginning  in  1971  was  due  to  major  receiver 
improvements,  including  the  installation  of  a  dual  IF  system. 

In  the  upper  halves  of  Fig.  2  we  used  black  curves  to  show  subjective 
representations  of  the  3-mm  data  (showing  the  actual  data  points  would  have 
overcrowded  the  figure);  we  used  dashes  where  the  data  are  sparse.  He  have 
been  cautious  in  drawing  these  representations  —  the  DR  21  results  indicate 
that  there  can  be  spurious  variations  of  2  -  3  Jy  amplitude  in  our  results 
(see  Sec.  II.C.lv).  So  the  rmader  may  make  an  unbiased  Interpretation  of  the 
data® ,  we  have  not  drawn  curves  through  the  actual  3-mm  data  in  the  lower 
halves  of  Fig.  2. 
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Table  I.  3.3-an  Observations8 


tKM  UK  CWTNK  TW  r«T  MO  HO  MTO.  TW  KM  fUff  OBBm.  MO  T*  n*TSTCM  STMDMO  BWK  onw  WM  «  N)  W8**vm 


16 


J 


IV.  CENTIMETER-WAVELENGTH  DATA 

These  3- am  results  ere  the  largest  available  set  of  aillimeter  monitoring 
data  with  high  temporal  resolution.  To  facilitate  comparisons  with  longer 
wavelength  data,  we  also  show  in  Fig.  2  the  four  sets  of  centimeter-wavelength 
data  which  have  the  most  nearly  comparable  resolution.  No  similar  body  of 
monitoring  data  exists  between  3  and  19  mm.  The  1.94-cm  data  are  from  Ref.  9 

and  private  communication  with  Dent;  the  2.8-  and  4.5-cm  data  are  from  Ref.  8; 

the  3.8-cm  data  are  from  Refs.  10  through  12  and  private  communication  with 

Aller.  Aller's  flux  densities  were  multiplied  by  1.044,  the  factor  derived  by 

Ledden*  to  bring  into  agreement  the  Michigan  and  Haystack  scales  for  these 
data.  Altschuler  and  Wardle,  who  used  the  NRAO  3-element  interferometer,  did 
not  observe  3C  273;  and  we  have  not  used  their  3C  84  results  because  they  did 
not  agree  with  near-simultaneous  Haystack  and  scaled  Michigan  3C  84  data.  The 
disagreement  is  probably  accounted  for  by  uncertainties  in  the  gain-elevation 
curve  for  the  interferometer  at  small  zenith  angles. 12,13 

The  centimeter  data  are  shown  in  the  upper  halves  of  Fig.  2.  To  avoid 
overcrowding  the  l-o  error  bars  are  omitted  when  they  are  less  than  twice  the 
height  of  the  symbol  (the  absence  of  an  error  bar  does  not  necessarily  mean 
the  error  is  as  small  as  or  smaller  than  the  size  of  the  symbol). 

Because  there  is  so  much  information  in  Fig.  2,  we  sought  alternate 
presentation  modes  to  provide  a  succinct  way  to  see  the  time  and  frequency 
relationships.  The  3-dimensional  plots  in  Fig.  3  are  one  result.  These  plots 
are  sequences  of  spectra.  The  identical  data  for  each  source  are  shown  from 
two  aspect  angles.  There  is  some  loss  of  information  because  the  light  curves 
had  to  be  smoothed  and  then  sampled  at  uniform  intervals  (0.02  years)  to 
prepare  the  plots.  But  they  do  give  a  "global”  impression  of  the  behavior  of 
the  sources.  Caution:  this  presentation  mode  does  not  indicate  the 
uncertainties  in  the  data. 
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Figure  3 


Three-diaenslonal  plots  of  the  data  froa  Fig.  2.  Curves  have  been 
drawn  by  eye  through  the  data;  the  curves  have  then  been  saapled  at 
Intervals  of  0.02  years  to  obtain  the  data  for  these  plots.  These 
three-dlaenslonal  presentations  aaount  to  sequences  of  spectra. 

The  saae  data  are  shown  froa  two  aspect  angles.  Caution:  these 
plots  do  not  portray  the  uncertainties  In  the  data. 


V.  DISCUSSION 


A.  POSSIBLE  CONNECTIONS  BETWEEN  EVENTS  AT  DIFFERENT  WAVELENGTHS 

Table  II  lists  by  date  and  wavelength  all  of  the  outbursts  (defined  as  a 
rise  and  fall)  that  we  have  discerned  in  Fig.  2.  The  table  also  Includes 
three  sharp  drops  (see  Sec.  V.E).  We  have  studied  Fig.  2  extensively  to 
discern  aore  than  just  the  obvious  possible  connections  between  events  at 
different  wavelengths.  Possible  connections  are  identified  by  entering  the 
dates  of  the  events  on  the  saae  line  in  Table  II.  The  criteria  used  in 
identifying  connections  included  time  coincidence  or  time  lag  (we  assumed  that 
the  conventional  adiabatic  expansion  model  pertains  and  thus  we  excluded  the 
possibility  of  events  at  longer  wavelengths  occurring  earlier),  similarity  of 
the  shapes  of  the  light  curves,  and  relative  amplitude.  Ue  considered  links 
only  between  the  events  nearest  in  time  —  i.e.,  we  did  not  consider  the 
possibility  of  “leapfrogged”  events. 

Caveats :  Subjectivity  is  Involved  in  identifying  events;  even  aore 
subjectivity  is  involved  In  making  connections.  This  subjectivity  emphasizes 
the  need  for  signatures,  such  as  Stokes  parameters,  in  addition  to  total 
intensity  as  a  function  of  time,  to  identify  the  same  outburst  at  different 
wavelengths.  Furthermore,  it  is  possible  that  two  closely-spaced  3-aa  events 
may  merge  at  longer  wavelengths  Into  a  single  event.  It  is  also  possible  that 
some  events  observed  at  3  am  may,  at  oentiaeter  wavelengths,  simply  have  such 
reduced  amplitude  that  they  are  unidentifiable. 

Perhaps  the  most  striking  feature  of  Fig.  2  is  that  the  3-am  and  centi¬ 
meter  light  curves  display  everything  from  high  correlation  to  almost  none. 
There  is  much  activity  in  the  3C  84  3-aa  data,  but  no  obvious  corresponding 
variability  at  centimeter  wavelengths.  The  mean  of  the  envelope  containing 
the  3-am  measurements  was  at  the  saae  flux  density  levels  as  the  centimeter 
data  in  the  late  1960s,  but  decreased  to  about  half  the  centimeter  signal 
level  by  1974,  For  3C  120  there  are  possible  correlations  between  some  events 
at  3  am  and  centimeter  wavelengths;  the  mean  of  the  envelope  of  the  3-am  data 
is  of  about  the  saae  intensity  as  the  oentiaeter  envelope  mean.  For  OJ  287 


21 


Table  II.  Dates  of  Peaks  of  Possibly-Connected 
Events  at  Different  Wavelengths* 

Wavelength  (an)  Remarks 


3.3  19.4  28  38  45 


—  1968.0: 

1968.6  —  1968.6?  —  Very  broad 

1969.1  1969.4? 

1969.65 

1969.4:  1969.6 ?:  1969.95  1970.1: 

1969.9:  1970.1?  1970.19? 

1970.3:  —  1970.5?  1970.65?:  1970.85?? 

1971.4  —  — 

*1971.85  1972.2?:  1972.2:  1972.2:  1972.3: 

1972.4  —  —  —  — 

1973.05:  — 

*!!«**?}  1973.6:  1973.75: 

1973.34  J 

1974.1 

3C  120 


1967.9 

1967.9 

1968.1 

1968.36: 

— 

— 

— 

*1968.67: 

1968.65: 

1968.66: 

1968.74: 

1968.85: 

— 

1969.11? 

— 

*1969.35: 

1969.6 

1969.6 

1969.6 

1969.7? 

1970.25: 

1970.3 

1970.3 

1970.4 

1970.75 

1970.8?: 

1970.9: 

1970.9 

1971.0 

1971.55: 

1971.6: 

1971.6: 

1971.7 

1971.6 

*1972.7 

1972.7: 

19728 

1972.85 

1972.9 

*1973.05: 

1973.2: 

1973.2: 

1973.2 

1973.2 

^Meaning  of  symbols: 

?  Existence  of  peak  uncertain. 

:  Date  uncertain  because  peek  very  broad  and/or  poorly  defined. 

—  Mo  reasonably  connected  event  found. 

*  Event  used  in  k- fact or  analysis  (Sec.  V.C). 

Mo  entry  In  the  table  Indicates  insufficient  or  no  data. 
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Table  II.  (Continued) 


Wavelength  (am) 


3.3 

19.4 

28 

38 

45 

1974.4: 

1974.4 

1974.4 

1975.3: 

1975.2: 

1975.25 

1975.55 

OJ  287 

1971.03: 

1971.03?: 

1971.2: 

1971.25: 

— 

1971.25?: 

1971.6?: 

1971.85 

— 

— 

*1972.12: 

1972.2: 

1972.2 

1972.2: 

1972.18 

1972.22? 

1972.3 

— 

1972.8 

— 

— 

— 

*1973.05: 

1973.1: 

1973.2 

1973.2: 

1973.5: 

1973.8 

1973.9 

1975.0: 

1975.0 

1975.4 

1975.44 

3C  273 

*1966.6 

1966.8: 

1966.8: 

<  1967.8: 

1968.2: 

1968.2: 

— 

— 

1969.3?: 

1969.5 

1969.6: 

1969.7: 

1969.55 

— 

1969.83 

1970.2? 

1970.3: 

— 

— 

— 

1970.7?: 

— 

— 

— 

1970.7?: 

*1971.8: 

1971.8: 

1971.8: 

— 

1972.0?: 

— 

— 

— 

1972.3: 

— 

1972.4: 

— 

— 

— 

— 

*1973.5: 

1973.6: 

1973.6 

1973.9: 

>1975.5 

Me 

1975.6: 

1975.8: 

BL  Lac 

_ 

1968.75: 

1968.75 

1968.8: 

1969.02 

1969.10: 

1969.12: 

1969.2: 

1969.2 

1969.19 

— 

— 

1969.48 

1969.49 

*1969.55: 

1969.65 

1969.66 

1969.67: 

*1969.9: 

1969.9: 

1969.96: 

1969.96 

*<1970.1: 

1970.2: 

1970.2: 

1970.2 

1970.21 

1970.35?: 

1970.35? 

1970.4?: 

1970.37: 

— 

1970.7? 

1970.7 

1970.65?: 

1970.67 

*>1970.86: 

aw 

1970.9?: 

1970.88 

<1970.95 

Me 

1970.9: 

Remarks 


Sharp  drop 


Sharp  drop 
Sharp  drop 

Very  broad 

Very  broad 
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Table  II.  (Continued) 


Wavelength  (na) 


Renarks 


3.3 

19.4 

28 

38 

45 

1971.1? 

1971.10: 

1971.08?: 

1971.10? 

— 

— 

1971.20 

— 

1971.35?: 

1971.35? 

*1971.45:: 

1971.5: 

1971.49 

1971.49 

1971.5 

*1971.83: 

<1971.9 

1971.83 

1971.9: 

1971.85: 

1972.1 

'"1972.15?: 

1972.14: 

<1972.15 

1972.3: 

1972.46? 

— 

~  — 

1972.6 

1972.62 

1972.6 

1972.60? 

— 

— 

1972.74?: 

1972.8 

1972.9 

1972. S3?: 

1972.93? 

1973.1 

1973.12 

— 

*1973.5 

1973.4: 

1973.52 

1973.55: 

1973.53 

1973.64 

1973.7? 

1973.72 

1973.7 

1973.9 

1973.9 

1974.4 

1974.42: 

1974.43 

— 

1974.63: 

— 

1975.23 

1975.25?: 

1975.43?: 

1975.42?: 

>1975.55  > 

1975.6 

1975.6 

24 


there  are  sone  correlations  between  3-mm  and  centimeter  events.  The  3-tun 
envelope  lies  above  the  envelope  of  the  centimeter  observations  in  1971  and 
1972,  but  is  at  the  level  of  the  centimeter  observations  in  1973  and  1974.  In 
3C  273  there  are  millimeter  to  centimeter  correlations  at  1966.6,  1971.8,  and 
1973.5;  there  Is  also  much  3-mm  activity  without  obvious  corresponding 
centimeter  activity.  The  mean  of  the  envelope  of  the  3-mm  measurements  has 
approximately  half  the  intensity  of  the  centimeter  envelope  mean.  For  BL  Lac 
there  is  clear  correlation  between  both  envelopes  and  events  at  3-mm  and 
centimeter  wavelengths.  For  all  five  sources,  generally,  the  flatter  the 
spectrum  from  3  to  45  mm,  the  closer  the  correlation  of  activity. 

For  six  of  the  BL  Lac  outbursts  the  spectral  index  was  ■  0  between  3  and 
45  mm.  Standard  models  would  suggest  that  the  optical  depth  from  3  to  45  mm 
was  very  small.  Events  were  thus  simultaneously  observable  at  all 
wavelengths.  Conversely,  for  those  objects  where  the  correlation  is  poor  or 
absent,  the  optical  depth  at  centimeter  wavelengths  may  be  high  enough  that  we 
don't  see  through  the  source  to  where  the  3-mm  activity  is  occurring.  In 
other  words,  the  3-mm  and  centimeter  emission  are  uncorrelated  because  they 
originate  from  different  regions  within  the  source. 

The  absence  of  correlation  between  the  3-mm  and  centimeter  intensity 
variations  in  many  Instances  is  in  contrast  to  the  high  correlation  between 

P 

short  and  long  centimeter  variations  found  by  Andrew  et  al.  and  the 
conclusion  by  Robbs  and  Dent  (1977)  of  high  correlation  between  3-mm  and 
longer  variations.  (The  approximately  85  sources  considered  by  Andrew  et  al. 
and  the  ten  considered  by  Hobbs  and  Dent  include  the  five  observed  for  this 
report.)  That  the  Hobbs  and  Dent  conclusion  is  based  on  undersampled  data  is 
clear  from  inspection  of  the  lower  half  of  Fig.  2,  where  their  3-mm  data  (red 
triangles)  are  plotted  with  our  3-mm  light  curves.  Their  data  are  simply  too 
sparse  to  catch  many  of  the  3-no  outbursts.  The  shortest  wavelength  data  that 
Andrew  et  al.  examined  is  the  19.4-mm  data  of  Dent  et  al. which  we  show  in 
the  upper  halves  of  Fig.  2.  Yet  we  often  do  not  find  any  correlation  between 
the  3-  and  19.4-mm  light  curves.  Thus  the  correlation  found  at  centimeter 
wavelengths  often  breaks  down  between  19  and  3  mm,  at  least  for  the  five 
sources  discussed  here. 
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B. 


COMMENTS  ON  SELECTED  EVENTS 


To  aid  the  reader,  we  draw  attention  to  the  following  interesting  events. 

3C  120.  Two  3-am  outbursts  peak  at  1972.2  and  1973.0.  The  associated 
centiaeter  outbursts  are  of  nearly  equal  amplitude,  while  the  3-am  outbursts 
have  greatly  different  aaplltudes. 

VLSI  data1*  indicate  a  super luainal  expansion  of  source  components  with  a 
nominal  aero  separation  tlae  of  1972.4  ±  0.5.  A  major  outburst  in  the  3 -mm 
and  centimeter  light  curves  began  at  1972.4;  this  3-mm  outburst  is  one  of  the 
two  or  three  strongest  3 -am  3C  120  events.  The  next  zero-separation  time 
shown  by  the  VLBI  data,  1974.4  ±  0.2,  coincides  with  a  small  centimeter 
outburst  and  a  possible  3-mm  outburst  (the  3-mm  data  are  limited). 

OJ  287.  The  strongest  OJ  287  3-mm  outburst  in  our  data  peaked  at  1972.17 
±  0.01;  strikingly,  this  outburst  was  terminated  by  one  of  the  sharpest  drops 
observed  at  3  am,  from  16  ±  2  to  7  ±  2  Jy  in  £  7  days  (see  Table  III).  A 
partial  recovery  occurred  during  the  next  several  weeks.  There  is  an 
indication  of  a  drop  at  28  mm  at  1972.22,  but  the  data  are  too  sparse  to 
discern  whether  the  decline  was  as  rapid  or  as  large  as  that  at  3  mm.  There 
are  at  least  three  oscillations  of  ~  20X  amplitude  at  38  nm  during  the  first 
three  months  of  1972  (see  the  insert  in  Pig.  2  of  Aller  and  Ledden23  for  an 
expanded-scale  plot  of  the  38-mm  data).  The  descent  portion  of  one  of  these 
oscillations  began  between  1972.16  “  29  February  and  1972.17  *  3  March,  l.e.,. 
at  about  the  same  time  as  the  drop  at  3  am.  No  further  38-mm  data  were 
obtained  until  15  days  later,  18  March,  at  which  time  the  signal  was  at 
another  oscillation  peak.  Because  the  *  502  drop  at  3  nm  was  extraordinary, 
whereas  the  observed  portion  of  the  decline  at  38  nm  was  not,  one  might 
speculate  that  the  38-am  signal  actually  declined  much  further  during  the  15- 
day  gap  in  observations. 

There  is  a  16Z  drop  in  the  28-mm  data  at  1972.07  “  26  January  (see  Table 
III).  The  expanded-scale  plot  of  these  data  (Fig.  3,  Ref.  17)  shows  that  the 
drop  occurred  within  1  day.  We  do  not  know  whether  the  decline  continued — it 
was  observed  on  the  last  day  of  a  multi-wavelength  monitoring  campaign.  The 
28-nm  measurements  did  not  resume  for  7  weeks.  Kikuchi  at  al.*®  observed  a 
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decrease  of  •  202  in  QJ  287'  s  72-mm  signal  over  100  min  on  8  February  1973 
(-1973.11)— see  Table  III. 


3C  273.  The  most  surprising  features  are  the  two  sharp  drops  at  3  bob: 
from  -  36  to  -  26  Jy  within  4  days  at  1969.55,  and  then  after  20  days, 
recovery  to  *>  38  Jy  within  6  days;  and  from  ■  41  to  -  20  Jy  within  2  days  at 
1969.83  (with  gradual  recovery  over  the  next  few  months) — see  Table  111. 

Three  observations  were  made  during  the  *  20-day  duration  of  the  first  drop; 
each  observation  consists  of  10  hours  of  integrations  spread  out  over  -  5 
day 8.  Because  the  telescope  was  in  operation  on  an  almost  round-the-clock 
basis,  it  is  very  unlikely  that  a  systematic  error  was  in  effect  for  20  days 
and  yet  had  no  observable  effect  on  any  of  the  several  other  objects 
monitored.  Scrutiny  of  the  original  data,  including  the  telescope  diagnostic 
information  recorded  every  25  min  (command  and  actual  coordinates, 
temperature,  relative  humidity,  refraction  computation  data,  antenna  tower 
tilt  data,  and  pointing  offsets),  uncovered  no  anomalies.  He  can  rule  out 
Sun-inrsldelobe  effects  because  the  drops  occurred  when  the  Sun  was  far  from 
3C  273  (they  were  within  15  deg  of  each  other  during  the  Interval  1969.71  - 
1969.79);  furthermore,  our  extensive  Mercury  observations  show  that  the  Sun- 
inrsldelobe  effects  for  our  antenna  do  not  maintain  the  same  amplitude  and 
sign  as  the  separation  angle  changes. 

There  is  no  clear  evidence  of  corresponding  drops  in  the  centimeter 
data.  For  clarity,  the  1969  and  1970  38-mm  data  are  shown  with  an  expanded 
intensity  scale  in  the  insert  in  the  lower  half  of  Fig.  2.  Only  Michigan  data 
are  shown  because  they  are  more  precise  and  more  frequent  than  the  Haystack 
38-mm  data;  when  error  bars  are  not  shown,  they  are  approximately  the  size  of 
the  symbols.  [Efanov _et_ ml.*®  found  declines  of  -  152  and  *  202  in  the  13.5- 
mm  flux  density  of  3C  273  in  only  -  4  hours  on  22  and  29  March  1976, 
respectively;  on  both  occasions  the  partial  recoveries  took  1  to  2  days. 
Gorshkov  et  al.2^  reported  a  3C  273  decrease  of  -202  at  35  mm  from  23  to  27 
May  1969  (1969.40),  but  the  published  observational  data  are  limited.] 

Linear  extrapolation  of  VLBI  data14  yields  1967.6  ±  0.4  as  the  start  time 
of  a  superluminal  expansion.  This  time  does  not  coincide  with  the  late  1965 
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time  of  the  strongest  3-mm  outburst  ever  observed  In  3C  273.  But  it  does 
coincide  with  a  3-nm  outburst  in  late  1967  (inferred  from  the  3-mm  decline 
observed  in  very  late  1967/early  1968)  and  a  centimeter  outburst  which 
commenced  in  late  1967.  This  3-mm  outburst,  though  poorly  defined  by  the 
data,  seems  to  be  the  second  strongest  event  observed  in  3C  273. 

BL  Lac.  The  coherence  of  the  light  curves  of  this  source  is  the  closest 
of  the  five  sources  presented  here.  The  time  lags  between  the  peaks  at  suc¬ 
cessively  longer  wavelengths  are  very  small  (days  to  weeks).  Although  this 
behavior  can  be  explained  by  the  simple  instantaneous-injection  expanding 
source  model,  more  complex  models  can  also  explain  it. 

C.  RELATIONSHIP  BETWEEN  OUTBURST  AMPLITUDE  AND  WAVELENGTH 

We  have  examined  the  data  of  Fig.  2  for  the  relationship 

ASX  a  Xk  (1) 

where  AS^  is  the  maximum  amplitude  above  the  quiescent  level  of  an  outburst  at 

wavelength  X.  For  a  simple  expanding  source  model  k  »  -1.0  (Ref.  25). 

o  +o  o 

However,  Andrew  et  al.  found  <k>  ■  -0.4  „q*15  for  44  sources  for  which  they 
had  28-  and  45-nn  data.  Andrew  et  al.  used  an  objective,  rote  technique  (see 
footnote  d,  Table  IV)  to  find  AS^  because  of  the  difficulty  in  choosing  a 
quiescent  level  (it  is  usually  assy  to  define  with  confidence  the  maximum  flux 
density);  their  rote  technique  is  probably  satisfactory  because  they  had  a 
fairly  large  sample.  Because  we  have  only  five  sources,  we  chose  to  make  more 
thorough  use  of  the  data  by  estimating  a  quiescent  level  for  each  outburst. 
However,  because  the  estimation  is  often  quite  arbitrary,  we  used  two 
different  procedures:  a)  linear  interpolation  between  the  troughs  on  both 
sides  of  the  outburst,  and  b)  extrapolation  forward  in  time  of  the  decay  curve 
of  any  preceding  outburst  and  backward  In  time  of  the  growth  curve  of  any 
succeeding  outburst.  The  outbursts  we  analysed  are  indicated  by  an  asterisk 
In  Table  II.  (These  practical  complications  make  it  difficult  to  compare 
observations  with  synthetic  light  curves  generated  by  models  unless  the 
synthetic  curves  Include  noise  and  overlapping  events.) 
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Table  IV.  <k>  Values8 


Procedure  used 

Wavelengths 

considered 

(mm) 

Interpolation** 

Extrapolation6 

Roted 

3.3  and  28 

-0.29  ±  0.33 

-0.31  ±  0.15 

-0.44  ±  0.15 

(n  -  18) 

(n  -  14) 

(n  -  5) 

3.3  and  38 

-0.35  i  0.24 

-0.37  ±  0.18 

-0.44  *  0.09 

(n  -  16) 

(n  -  12) 

4A 

1 

0 

3.3  and  45 

-0.39  ±  0.23 

-0.44  *  0.18 

-0.45  *  0.03 

(n  -  15) 

(n  -  13) 

(n  -  5) 

8  The  statistical  scatter  Is  Indicated  by  the  cited  standard  deviations.  The 
number  of  connected  events  analyzed,  n,  is  given  in  parentheses. 

b  Outburst  quiescent  levels  determined  by  linear  interpolation  between 
adjoining  troughs  In  the  light  curves. 

c  Outburst  quiescent  levels  determined  by  extrapolation  of  the  decay  curves 
of  any  preceding  outbursts  and  of  the  growth  curves  of  any  succeeding  out-: 
bursts.  Fewer  events  are  suitable  for  quiescent  level  determination  by 
this  procedure  than  by  the  linear  interpolation  procedure. 

d  The  rote  procedure  of  Andrew  _et^  al .  (1978)  wherein  individual  outbursts  are 
not  considered.  The  peak  and  quiescent  levels  are  simply  set  equal  to  the 
maximum  and  minimum  flux  densities,  respectively,  observed  over  the  entire 
monitoring  period.  Justification  for  this  procedure  is  given  by  Andrew  et 
el. 
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He  calculated  k  for  three  wavelength  Intervals:  from  3.3  to  28,  38,  and 
45  mm.  The  values  of  <k>  for  each  Interval  and  for  the  two  different  quies¬ 
cent  level  determination  procedures  are  contained  in  Table  IV.  The  two  proce¬ 
dures  give  essentially  the  same  <k>  values,  although  the  values  determined  by 
the  potentially  more  realistic  extrapolation  procedure  have  smaller  standard 
deviations.  The  fourth  column  gives  the  values  determined  by  the  Andrew  et 
al.  rote  procedure;  they  are  not  significantly  different  from  the  other 
values.  We  do  not  understand  the  surprisingly  small  scatter  In  the  rote  pro¬ 
cedure  results;  perhaps  it  is  fortuitous.  To  obtain  an  indication  of  the  sen¬ 
sitivity  to  the  adopted  quiescent  level,  we  computed  k  values  with  quiescent 
levels  differing  by  ±  52;  the  resulting  changes  in  the  <k>  values  were  smaller 
than  the  uncertainties  due  to  statistical  scatter.  From  the  Table  we  see  that 
for  the  intervals  3.3  -  28,  3.3  -  38,  and  3.3  -  45  on,  <k>  »  -0.4  (std.  dev.  * 
0.2),  in  excellent  agreement  with  the  Andrew  et  al.  value  for  44  sources  for 
28  -  45  mm.  Had  our  measurements  of  the  3-mm  outburst  amplitudes  been,  on 
average,  larger  or  smaller  by  a  factor  of  two  (a  very  large  error),  we  would 
have  obtained  <k>  values  of  -0.7  or  -0.1,  respectively. 

Caveats.  Random  errors  may  be  present  if  we  have  made  wrong  connections 
between  3-mm  and  centimeter  outbursts.  Some  3-mm  events  may  have  been  left 
out  of  the  analysis  because  we  were  unable  to  identify  their  corresponding 
centimeter  events;  thus,  we  may  be  inherently  biased  toward  less  negative 
values  of  k. 

Andrew  et  al.  also  analysed  23  sources  for  which  Altschuler  and  Wardle1^ 

had  37-  and  111-an  data  and  found  <k>  ■  -0.4  ±  0.15.  Thus  we  have  the 

-0  4 

remarkable  result  that  outburst  amplitude  varies  approximately  as  X  *  from 
3.3  to  111  mm,  a  wavelength  ratio  of  34.  [This  ratio  is  actually  slightly 
larger  because  the  emission  observed  at  3.3  ma  is  emitted  at  the  source  at 
3. 3/(1  +  z)  am.  We  have  neglected  this  relativistic  correction  because  the 
largest  s-value  of  our  five  sources  is  only  0.3  (0J  287).]  We  know  of  no 
theory  which  predicts  such  a  weak  dependence  of  outburst  amplitude  on  wave¬ 
length.  Therefore  this  result  nay  place  a  strong  constraint  on  models  of  the 
variable  emission  of  extragalactlc  radio  sources. 
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TIME  LAGS  BETWEEN  OUTBURSTS  AT  DIFFERENT  WAVELENGTHS 


On  the  assumption  that  the  possible  connections  in  Table  II  ate  real,  we 
have  calculated  the  tine  lags,  if  any,  between  outbursts  at  3  am  and  longer 
wavelengths.  Although  the  epochs  of  event  peaks  are  often  poorly  determined , 
we  have  treated  then  as  accurate  and  precise,  even  when  given  to  only  one 
decimal  place  and  when  uncertainty  is  indicsted,  sinply  because  no  better  data 
axe  available.  Figure  4  shows  the  average  tine  lag  between  possibly-connected 
events  as  a  function  of  the  difference  between  observing  wavelengths.  The 
vertical  bars  indicate  the  dispersion  in  values,  not  the  standard  error  of  the 
aean;  the  nuaber  of  events  in  the  saaple  is  also  indicated.  The  time  lags  are 
smallest  in  BL  Lac  (~  0.06  y)  and  longest  in  3C  84  (~  0.2  to  ~  0.6  y).  Within 
the  limitations  of  the  data,  there  is  a  suggestion  that  longer  tine  lags  are 
associated  with  larger  wavelength  Intervals.  Note  that  the  two  Seyfert 
galaxies  in  our  sample  (3C  84  •  NGC  1275  and  3C  120)  have  longer  time  lags 
than  the  two  Lacertids  (OJ  287  and  BL  Lac). 

E.  SHARP  DROPS  (QUENCHINGS) 

We  have  found  in  the  Aerospace  3-arn  data  no  outbursts  on  time  scales  of 
one  to  a  few  days.  Neither  did  we  find  such  rapid  3 -am  outbursts  in  several 
searches  with  the  NRA0  11-m  antenna  for  Intraday  and  lnterday  variations  of 
extragslactic  sources^.  But  we  have  found  in  the  Aerospace  and  NRA0  data  a 
nuaber  of  rapid,  large  amplitude  (>  302)  drops,  or  quenchings,  on  time  scales  _ 
of  one  to  several  days.  Table  III  contains  information  on  quenchings  at 
various  radio  wavelengths  In  these  and  other  sources;  the  tabulation  is  not 
meant  to  be  exhaustive.  Three  quenchings  hsve  time  scales  of  just  a  few 
hours. 

Caveats .  The  definition  of  quenching  and  the  selection  of  quenching 
events  is  arbitrary.  We  selected  as  quenchings  those  drops  which  are  far  more 
rapid  than  other  increases  and  decreases  observed  In  the  given  source  at  the 
given  wavelength.  We  did  not  Include,  for  example,  any  of  the  three  0J 
287  drops  of  ~20%  at  38  ■  shown  In  the  Insert  of  Fig.  2  of  Baf.  23  because 
their  data  show  that  ~  *202  variations  were  common  at  that  time  in  0J  287. 

Had  there  been  only  one  Isolated  rapid  drop  we  probably  would  have  included 
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Figure  4  The  avenge  time  lag  between  possibly-connected  events  as  a 

function  of  the  difference  between  observing  wavelengths.  The 
vertical  bars  indicate  the  dispersion  in  values,  not  the  standard 
error  of  the  mean;  the  numerals  indicate  the  number  of  events  in 
the  samples. 
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it.  To  the  extent  that  similar  first  derivatives  may  indicate  similar 
physics,  our  selection  criterion  may  be  causing  us  to  omit  some  physically 
interesting  events.  The  sharpness  of  these  drops  is  in  narked  contrast  to  the 
relatively  slow  decay  of  outbursts  predicted  by  all  published  models  we  know 
of.  The  observations  of  quenchings  thus  place  a  new,  and  possibly  strong, 
constraint  on  models  of  extragalactic  sources.  Perhaps  eclipsing, 
obscuration,  or  rotation  phenomena  are  involved. 
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VI.  SUMMARY 


Mine  years  of  3-sn  observations  of  variable  extra  galactic  radio  sources 
have  been  presented  and  compared  with  centimeter-wave  observations. 

1)  Useful  information  on  extragalactlc  variable  sources  can  be  obtained 
with  an  antenna  of  only  modest  size. 

2)  Significant  information  on  variations  of  extragalactlc  radio  sources 
at  short  millimeter  wavelengths  is  lost  if  temporal  resolution  of  about  one  to 
three  days  is  not  achieved. 

3)  The  degree  of  correlation  between  3-an  and  centimeter-wave  variations 
ranges  from  high  to  low  or  almost  non-existent.  In  general,  the  flatter  the 
spectrum,  the  better  the  correlation. 

4)  The  relationship  between  time  lag  and  wavelength  for  possibly- 
connected  outbursts  appears  to  be  characteristically  different  for  different 
objects. 

5)  No  pronounced  3-mm  outbursts  on  a  time  scale  of  one  to  a  few  days  were 
observed,  but  three  sharp  3-an  drops  or  quenchings  were.  These  Aerospace 
measurements,  plus  our  measurements  with  the  NRAO  11-m  antenna,  suggest  that 
rapid  drops  at  3  am  may  be  more  common  than  rapid  rises. 

k 

6)  The  amplitudes  of  outbursts  vary  as  X  ,  where  <k>  ■  -0.4,  from  3.3  to 
111  am,  a  factor  of  34. 

Results  5  and  6  may  place  important  new  constraints  on  physical  models  of 
extragalactlc  variable  radio  sources. 
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